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We investigate the superconducting (SC) state in two-dimensional extended Hubbard 
model at a quarter-filling with the on-site repulsive interaction (U) and the nearest- 
neighbor one (V), which lead to the spin fluctuation and the charge fluctuation, respec- 
tively. The effect of the charge fluctuation on the onset of the SC state is calculated 
by estimating the pairing interaction within the random phase approximation. Solving 
the linearized Eliashberg equation with the anomalous self-energy as the function of 
both momentum and frequency, we examined the singlet SC state on the plane of U 
and V with the fixed temperature. We found a novel result that the SC state with the 
d xy symmetry retains for most parameters of U and V, although the pairing interac- 
tion changes the sign from repulsive one into attractive one with increasing V. The 
reentrant transition as the function of temperature occurs close to the onset of the SC 
state. 
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Organic conductors exhibit the variety of superconducting (SC) states, as found in quasi- 
one-dimensional (Q1D) materials of Bechgaard salts, 1 ) or quasi-two-dimensional (Q2D) ma- 
terials of ET salts. 2 ' 3 -* It has been often asserted that the spin fluctuation is the origin of 
the SC state since the spin density wave state is located next to the SC state in the phase 
diagram. However there are some ET salts suggesting another possibility that the charge 
fluctuation could be the origin of the SC state. In #-ET 2 X salts, where the dihedra angle 
between two adjacent molecules changes with the choices of the anion X, it has been shown 
that the state with the charge ordering moves into the SC state with decreasing the dihedral 
angle 2 -* corresponding to the increase of transfer energy. It is noteworthy that this SC state 
is located besides the charge ordering state. Moreover there is the recent observation of the 
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SC state in the Q2D conductor, f3' -ET 2 IC1 2 , which appears under the pressure of ~ 8 GPa. 4 ^ 
This indicates the increase of the nearest-neighbor repulsive interaction in addition to the 
increase of the transfer energy. 

In 2D systems, there are several theoretical work for the SC state, which is based on 
the FLEX ( fluctuation exchange ) approximation. 5 ) Especially, the spin fluctuation as the 
origin of the singlet SC state has been calculated for several transfer energies. By using a 
model with the on-site repulsive interaction U and a variation from the square lattice to the 
triangular lattice, 6 ) it has been shown that the SC state at half-filling is mainly determined 
by the state with the d x 2_ y 2 symmetry. On the other hand, in the presence of the nearest 
neighbor repulsive interaction V, there appears the SC state with the d xy symmetry due to 
the charge fluctuation, as shown by the RPA (random phase approximation) 7 ) with a choice 
of U and V and by the slave-boson theory with a quarter-filled band. 8 ) However, it is not 
yet clear how the interplay of U and V leads to a crossover from the spin fluctuation to the 
charge fluctuation especially for the quarter-filled band. 

In the present paper, such a role of charge fluctuation is studied by using the two- 
dimensional extended Hubbard model at the quarter-filling of hole doping. The case of 
the square lattice is examined as the first step since the transfer integral in organic conduc- 
tors is very complicated. Within RPA, 9 ' 10,7 ) we examine a phase diagram of SC, SDW (spin 
density wave) and CDW (charge density wave) on the plane of U and V and also these states 
as the function of temperature. 

We consider the expended Hubbard model given by 



where quantities U and V are the coupling constants for on-site and nearest-neighbor re- 
pulsive interactions. The quantity denotes the band energy of the 2D square lattice, 
£k = —2t(cos k x a + cos k y a) — /i, where t is the nearest-neighbor transfer energy and fi is 
determined to obtain the 3/4-filled band, i.e., 1/4-filled hole band. In eq. (1), Nl is the total 
number of the lattice site and a is the lattice constant. We set t — 1, a — 1 and &b = 1- 
First, we calculate spin and charge susceptibilities given by 



with T and T T being the temperature and the time ordering operator of the imaginary time. 
The spin susceptibility is given by A = s and f a = 1(= —1) for a =| (=|) while the charge 
susceptibility is given by A = c and f a = 1. The irreducible susceptibility x°(q, w;), which 
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denotes eq.(2) without the interactions, is given by 

X (q 1 ^) = -FEG (k + q ) e m + Wi )G (k,g , (3) 

where io\ (e m ) is the Matsubara frequency for boson (fermion) and G°(k, e m ) = [ie m — £k] _1 - 
Treating U and V terms of eq. (1) within RPA, eq.(2) is calculated as 

X s (q,^) = X °(q,^)/(1 " U X °(<l,wi)) , (4) 

X c (q, WI ) = x°(q, w,)/(l + (C/ + 2V(q)) X °(q, w,)) , (5) 

where V(q) = 2l / (cosg x a + cos q y a). 

Next, we derive the pairing interaction, P(q, a;j), corresponding to the vertex for the singlet 
state 

P(q,a;,) / drc| c _ q , T (r)c t _ k+qa (r)c_ ki (0)c kT (0)e^ T , (6) 

which is used in the Eliashberg equation for the SC state. In terms of eqs. (4) and (5), the 
pairing interaction P(q, is obtained as 

P(q,^)=f/ + V(q) + ^V(q,^) 

_( * u* + 2 [/y(q) + 2\/(q) 2 ) X c (q, w,) . (7) 
The SC state is determined by the anomalous self-energy, S a (k, e m ), defined by 
£ |dr(T T c k , T (r)c_ k , i (0))e i ^-P(k-k / ,e m -e m ,) . 

(8) 



The linearized Eliashberg equation, which gives the onset of the finite value of E a (k, e m ), 
is expressed as 11 ) 

AS a (k,e m )= ]T ir(k,e m ;k / ,e m 0S a (k / ,e m , (9) 

k',e m / 

X(k, e m ; k', e m = -^P(k - k', e m - e m ,)\G°(k', e m ,)\ 2 . 

The SC transition temperature T c is given by the condition, A = 1. Since the characteristic 
momentum dependence of P(q, uJi) is essentially independent of u>i, and P(q, uii) decreases 
monotonously with \u>i\, eq. (9) is calculated by assuming a form given by 

£ a (k, em ) = £?(e m )£^(k) , (10) 

where E°(e m ) is evaluated by a sample average with some choices of k. We solve the linearized 
Eliashberg equation numerically by using an iterative method where the symmetry in k-space 
is rigorously determined without the assumption for the pairing symmetry. 
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Fig. 1. Spin (x s ) and charge (x c ) susceptibilities for (U, V) = (3.7, 0), (0, 0.8) where T = 0.01. 



We examine the role of the spin and charge fluctuations by varying U(> 0) and V(> 0). 
From eqs. (4) and (5), the spin and charge susceptibilities are calculated. For U ^ and 
V = 0, the spin fluctuation is much larger than the charge fluctuation as seen from the 
solid line in Fig. 1. This is because the magnetic state is expected by the repulsive on site 
interaction leading to the enhancement of the spin fluctuation, \ s ■ The nesting condition 
at 3/4-filling gives rise to \ s being dominant in the interval region of (q x ,Qy) = (^j 0) and 
(tt, ±7r), while the peak of \ s a t 1/2-filling appears at (tt,tt). For V U, the charge 
susceptibility, \ C i becomes much larger than \ s , since V induces CDW (or charge ordering) 
and then the enhancement of the charge fluctuation. In this shown by the dashed 

line in Fig. 1, the peak of x c appears as the ridge connecting two points of (tt, tt — b) and 
(tt — b, 7r) with b ~ it/2 in a quarter plane of (q x , q y ). This comes from the effect of V(q) in 
eq.(5) exhibiting a peak at (n,ir) in addition to that of the nesting condition. 

Based on Fig. 1, we calculate the pairing interaction, P(q, wi), which is given by eq.(7). 
For U 7^ and V = 0, P(q, 0) is always positive and shows the large magnitude around 
(tt, 0) with the peak located between (tt, 0) and (tt, tt) as seen from the solid line of Fig. 2. 
The q dependence of P(q, 0) is similar to that of X s (qj0) m Fig- 1 since eq.(7) leads to 
P(q, 0) oc x s (q> 0) for V = and large \ s ■ However, for V ^ 0, -P(q) becomes different from 
charge and/or spin susceptibilities since the second line of eq.(7) leads to P(q) having the 
sign opposite to x c - With increasing V/U , P(q) begins to decrease especially in the region 
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Fig. 2. Pairing interaction P(q,0) at T = 0.01 for (U, V) = (3.7,0) (solid curve), (3.0,1.15) (short-dashed 
curve), (2.0,1.1) (dot-dashed curve), (1.0,1.0) (dashed curve) and (0,0.8) (dotted curve). 



of (tt, 0) and (tt, tt) where the peak of P(q) is reduced noticeably and the dip around (tt, tt) 
grows leading to the attractive interaction. Thus it turns out that, with increasing V, the 
main part of pairing interaction becomes attractive and takes a broad dip around (tt, tt) but 
still remains repulsive in the region being away from (tt, tt). Apparently, such an attractive 
interaction indicates the s-wave pairing but our obtained result differs from such a pairing 
as shown later. 

The onset for the SC state is obtained by calculating eq.(9) numerically. The parameters 
U and V are calculated so that the maximum eigen value of A may become unity with the 
fixed T. In Fig. 3, the solid curve corresponding to the onset of the SC state is calculated 
for T = 0.01. On the dashed lines, either x s (cL)0) or X c (ci>0) diverges at T = 0.01, i.e., 
either the incommensurate SDW or the incommensurate CDW is obtained for (U, V) in the 
right hand side or the upper side of the respective line. The region enclosed by the solid 
line, the horizontal axis and the vertical axis corresponds to the normal state. The SC state 
is obtained in the region enclosed by solid line and dashed two lines. The critical value 
of V for the onset of the SC state increases with increasing f/(<3). This comes from the 
competition between the spin fluctuation and the charge fluctuation where the sign of the 
respective pairing interaction is opposite each other as shown in Fig. 2. Thus it is found that 
the pairing interaction corresponding to the onset of the SC state is mainly determined by 
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Fig. 3. Phase diagram on the plane of U and V at T — 0.01 where the superconducting state (SC) is 
obtained in the region between the solid curve and the dashed curves. The solid curve denotes the onset 
for the d-wave (d xy ) SC state and the dotted line denotes that for the s-wave SC state. In the inset, the 
detail close to U — is shown where solutions for both s (dotted line) and d xy (solid line) exist for U < 0.2. 



the charge fluctuation for U ^3.5 while the spin fluctuation is dominant for 3.5 %U ^$3.7, 
i.e., close to the upper bound of U. The SC region is enlarged with decreasing temperatures 
since the solid curve moves toward the origin (0,0) rapidly compared with that of the dashed 
line. 

In order to examine the symmetry of the order parameter of the SC state, k dependence 
of the quantity S a (k, mT), is shown in Fig. 4 for three kinds of parameters. The solid 
line, dashed line and dotted line correspond to (U,V) = (3.7,0), (1.0,1.0) and (0,0.8), 
respectively. In the case of (U, V) = (3.7, 0) and (1.0, 1.0), we found following two features 
: (1) the sign of S a (k, mT) in the region of k y > are different from that of k y < 0, (2) 
S a (k, mT) = on k x and k y axis and on the boundary of the Brillouin zone, indicating 
the relation, Ti a [k x ,ky,e m ) = —Yi a [k x ,—ky,e m ) and Ti a [k x ,ky,e m ) = — S a (— k x , k y , e m ). By 
noting the symmetry of the square lattice, this state corresponds to the d^-like symmetry. 
In the case of (0,0.8), on the other hand, the sign of S a (k, mT) remains unchanged except 
for the region near (jr, ir) and S a (k, mT) has a maximum near (0,0). Thus the dotted line 
of Fig. 4 indicates the s-like symmetry. In Fig. 3, such a s state is shown by the dotted line. 
The d xy state in the case of (U,V) = (3.7,0) is reasonable due to the repulsive interaction 
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located around q = (tt, 0) and seems to be consistent with that of Kondo 12 \ who treated 
the ?7-term perturbationally. For V being large with the moderate strength of U (dashed 
curve of Fig. 4), the SC state is still given by d xy although the paring interaction of Fig. 2 
is almost attractive. This is understood as follows. Both d xy state and s state can gain the 
the energy from the attractive interaction close to (tt, tt) although the latter is slightly larger 
than the former due to the absence of the node. However d xy state can be chosen due to the 
energy gain from the the repulsive interaction around (tt, 0), which remains for the case of 
the dashed curve. In the case of (U, V) = (0, 0.8), Fig. 2 shows that the repulsive interaction 
around (tt,0) becomes sufficiently small leading to s state.. The inset of Fig. 3 shows that 
d xy state still exists as a metastable state at U = 0. It is expected that the difference of the 
Helmholtz free energy between s state and d xy state is small for the case of small U. We 
note that there is no node for both d xy state and s state on the Fermi surface of the present 
case of the 3/4- filled band. 
Here, we analyze S a (k, mT) by using the Fourier transformation given by 

S a (k,i7TT)=^C mn mn (k) , (11) 

mn 

with <f) mn (k) = e mfea e m ^7(27r). In the case of (U, V) = (3.7,0), the ratio of the d xy - 
component (i.e., (fi dxy (k) = 2 sin k x a sin k y a) to other components with larger momentum is 
given by 1 : 0.50, while s- and cL^^-component are absent. In the case of (U, V) = (1.0, 1.0), 
the corresponding ratio is 1 : 0.37. For case of U ~ 0, such a ratio remains approximately 
the same, and the s- and d^^-components also do not exist, whereas the s-state is located 
very close to that of the the G^-state in the U-V plane. For the case of (U, V) = (0,0.8), 
the extended s state shows three components consisting of the s-component (0 s (k) = 1), 
an extended s-component (0 cxt _ s (k) = cos k x a + cos k y a) and other components with larger 
momentum, where the respective ratio is given by 1 : 0.59 : 0.39. We note that S a (k, mT) 
becomes negative near (vr,7r). 

Finally we examine the SC by varying temperature in Fig. 5. The frequency dependence 
in eq.(10) is discarded since such an effect becomes small with increasing temperature. With 
increasing T the critical values of V for onset of both the SC state and the CDW state 
increase at low temperatures (T^O.l) and at high temperatures (0.6 ^T), while the value 
of V decreases at the intermediate temperatures (0.1 ^ T < 0.6). Such a reentrant transition 
of the CDW state and the SC state as the function of temperature, which resembles that 
of reference 8, originates in the property of eq. (5). At low temperature, the maximum of 
eq.(5) is given by the nesting wave vector of x°(q, u;;), which is away from (tt, 7r) as shown 
by the dotted curve of Fig. 1. With increasing temperature, x°(q, u>i) as the function of q 
becomes broadened. The nesting wave vector begins to move toward (tt, tt) at T ~ 0.08 and 
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Fig. 4. Self-energy for the SC state , S(fc,i7rT) for (U,V) = (3.7,0) (solid line), (1.0,1.0) (dashed line) and 
(0,0.8) (dotted line) where T = 0.01. The inset shows the behavior in the quarter plane indicating that 
the sharp peak around (7r/2, 7r/2) for the d xy state and the broad peak around (0,0) for s state. 



also at the intermediate temperatures due to a competition of x°(q, w/) and |V(q)|, which 
results in the increase of eq.(5). At high temperatures, eq.(5) with a maximum at q = (tt, it) 
decreases due to the further broadening of x°(q, cj/). 

In summary, we have examined the effect of nearest neighbor repulsive interaction on the 
SC state for the 2D extended Hubbard model with a square lattice at 3/4- filling. We found 
the SC state with d xy symmetry for both cases of U and V, although the pairing interactions 
are repulsive for U and attractive for V. The spin fluctuation induced by U gives the SC 
state with the symmetry being the the same as that of the charge fluctuation induced by V. 
However the onset of the SC state needs larger magnitudes of interactions of U and V since 
these two fluctuations are incompatible as seen from Fig. 3. Thus it is interesting to extend 
the present study to the case of the SC state of organic conductors, in which exotic states 
may be expected due to the variety of transfer energies. 
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Fig. 5. Phase diagram on the plane of T and V, where U = 3. 
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